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Primordial magnetic fields will generate non-Gaussian signals in the cosmic microwave background (CMB) 
as magnetic stresses and the temperature anisotropy they induce depend quadratically on the magnetic field. We 
compute a new measure of magnetic non-Gaussianity, the CMB trispectrum, on large angular scales, sourced via 
the Sachs-Wolfe effect. The trispectra induced by magnetic energy density and by magnetic scalar anisotropic 
stress are found to have typical magnitudes of approximately a few times 10^^'' and 10~^^, respectively. Obser- 
vational limits on CMB non-Gaussianity from WMAP data allow us to conservatively set upper limits of a nG, 
and plausibly sub-nG, on the present value of the primordial cosmic magnetic field. This represents the tightest 
limit so far on the strength of primordial magnetic fields, on Mpc scales, and is better than limits from the CMB 
bispectrum and all modes in the CMB power spectrum. Thus, the CMB trispectrum is a new and more sensitive 
probe of primordial magnetic fields on large scales. 



Magnetic fields are ubiquitous in the Universe from planets 
and stars to galaxies and galaxy clusters yet the origin 

and evolution of large-scale magnetic fields remains a puzzle. 
A popular paradigm is that magnetic fields in collapsed struc- 
tures could arise from dynamo amplification of seed magnetic 
fields [2]. The seed field could in turn be generated in as- 
trophysical batteries 13i] or due to processes in the early uni- 
verse [stl- Indeed recent 7-ray observations claim to find 
a lower limit to an all-pervasive intergalactic magnetic field 
that fills most of the cosmic volume |6], which would per- 
haps favor a primordial origin. A primordial magnetic field 
can be generated at inflation fl, or arise out of other phase 
transitions in the early Universe (sj. As yet there is no com- 
pelling mechanism which produces strong coherent primor- 
dial fields. Equally, the dynamo paradigm is not without its 
own challenges in producing sufficiently coherent fields and 
sufficiently rapidly ||2(]. Therefore, it is useful to keep open the 
possibility that primordial magnetic fields originating in the 
early universe play a crucial role in explaining the observed 
cosmic magnetism. 

In this context it is important to investigate every observ- 
able signature of the putative primordial magnetic fields. Con- 
straints on large-scale primordial magnetic fields have already 
been derived using the cosmic microwave background (CMB) 
power spectrum IlZIIsI] and Faraday rotation [9]. However, the 
effects of a magnetic field on the CMB are relatively more 
prominent in its non-Gaussian coiTelations. This is because 
magnetic fields induce non-Gaussian signals at lowest order 
as the magnetic energy density and stress are quadratic in 
the field. On the other hand, the standard inflationary per- 
turbations, dominated by their linear component, can source 
non-Gaussian correlations only with higher-order perturba- 
tions and thus necessarily produce a small amplitude of CMB 
non-Gaussianity (cf. iiciil ill ). Primordial magnetic fields can 
induce appreciable CMB non-Gaussianity when considering 
the bispectrum fl2. 13]. Our previous calculation of the mag- 
netic CMB bispectrum sourced by scalar anisotropic stress led 
to a 2 nG upper limit on the primordial magnetic field's am- 



plitude on Mpc scales [Il4l]. However, higher-order measures 
of non-Gaussianity remain unexplored and as we show here, 
could be very useful to set further constraints on primordial 
magnetic fields. 

In this letter, we present the first calculation of the contri- 
bution to the CMB trispectrum induced by a primordial mag- 
netic field. In particular, we consider the magnetically induced 
Sachs-Wolfe effect sourced by a stochastic primordial mag- 
netic field. We show that the trispectrum does significantly 
better than the bispectrum in constraining the large-scale mag- 
netic field via CMB non-Gaussianity, considering both mag- 
netic energy density and magnetic scalar anisotropic stress as 
sources. This reveals a new and effective probe to investigate 
primordial magnetic fields on large scales. 

We consider a Gaussian random stochastic magnetic field 
B characterized and completely specified by its power spec- 
trum M{k). We further assume the magnetic field to be non- 
helical. On galactic and larger scales, any velocity induced 
by Lorentz forces is generally too small to appreciably dis- 
tort the initial magnetic field |15]. Hence, the magnetic field 
simply redshifts away as B(x, i) = 6o(x)/a^, where, 60 is 
the magnetic field at the present epoch (i.e. at 2 = or 
a — 1). We define b{k) as the Fourier transform of the mag- 
netic field bQ{x). The magnetic power spectrum is defined 
by the relation {b,{k)b*{q)) = {2TrfS{k-q)P,j{k)M{k), 
where Pij{k) — {6ij — kikj/k"^) is the projection operator 
ensuring V • bo = 0. This leads to (bg) = 2 / {dk/k)Al{k), 
where A^(A:) = k^M{k) / (27r^) is the power per logarithmic 
interval in k space present in the stochastic magnetic field. We 
assume a power-law magnetic power spectrum, M (fc) = Afc" 
that has a cutoff at A: = kc, where kc is the Alfven-wave damp- 
ing length scale Iil5i1 . We fix A by setting the variance of the 
magnetic field to be Bq, smoothed using a sharp /c-space fil- 
ter, over a "galactic" scale kc = Ih Mpc^^. This gives (for 
n > —3 and for k < kc) 
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The spectral index n is restricted to values close to and above - 
3, i.e., an inflation-generated field, as causal generation mech- 
anisms can only produce much bluer spectra ||l6ll. Further, 
blue spectral indices are strongly disfavored by many obser- 
vations like the CMB power spectra |7J. We choose to split 
the contributions to the CMB trispectrum into that sourced by 
magnetic energy density il b and by scalar anisotropic stress 
Hb rather than the com pen sated and passive magnetic pertur- 
bation modes of Ref. OlTll . The subdominant compensated 
mode is a linear combination ofils and Hb whereas the pas- 
sive mode is the Hb perturbation considered here. 

The Sachs-Wolfe type of contribution to the CMB temper- 
ature anisotropy induced by the energy density of magnetic 
fields Ii18h20i] . can be expressed as 
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Here, ^Bix) = B\x,t)/{87rp^{t)) = bl{x)/{87rpo), 
where p-y{t) and po are, respectively, the CMB energy den- 
sities at a time t and at the present epoch. In the same manner 
as the usual Sachs-Wolfe effect, the AT/T given above is for 
large angular scales. For numerical estimates we use the most 
recent estimate of Bonvin and Caprini (Eq. 6.12 of L2flJ) ex- 
pressed according to our definitions as TZ — — 0.2i?^/3 ^ 



-0.04 where 



0.6 is the fractional contribution of ra- 



diation energy density towards the total energy density of the 
relativistic component. The unit vector n is along the direc- 
tion of observation from the observer at position a;o and D* 
is the (comoving angular diameter) distance to the surface of 
last scattering. We have assumed instantaneous recombina- 
tion which is a good approximation for large angular scales. 

The temperature fluctuations of the CMB can be ex- 
panded in terms of spherical harmonics to give AT{n)/T = 
Y.im ahnYhn{n), where 



47r 



aw. 



_fk_ 

(2^ 



nnB{k)j,{kD*)Y;„M 0) 



Here, ilB{k) is the Fourier transform of ^^^(a;). Since 
fis(cc) is quadratic in bo{x), we have a convolution flsik) = 
(l/(27r)3) / d^s b,{k + s)b* {s)/ {St: pq). The trispectrum 



j,m-^mjn^m^ , or the four-point correlation function of the 
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CMB temperature anisotropy in harmonic space, in terms of 
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Eq. OJ we can express T ^ ^ 
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withCi234 - {nBik,)nBik,)nB{k,)nBikJ}. The four- 
point correlation function of ftB{k) involves an eight-point 
correlation function of the fields. Using Wick's theorem, for 
Gaussian magnetic fields, we can express the magnetic eight- 
point correlation as a sum of 105 terms involving the mag- 
netic two-point correlation. Neglecting the 45 terms propor- 
tional to 6{k) that vanish and the 12 terms proportional to 



6{ki + kj) that represent the unconnected part of the four- 
point correlation, we are left with 48 terms. A long calcula- 
tion involving the relevant projection operators gives C1234 = 
S{ki + fc2 + fcs + ^4) V-'i234' where ^^^3^ a mode-coupling 
integral over a variable s and also involves angular terms. The 
full expression for ^^1234 will be presented in our detailed pa- 
per |21]. For simplicity we evaluate the mode-coupling in- 
tegral ^/'i,34 in two cases: (I) considering only s-independent 
angular terms for all equal-sided configurations and (II) taking 
all angular terms for the collinear configuration. Considering 
s-independent terms only for a general configuration, we find 
"01234 = -8/(87r/9o)'' 2:where 
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We perform the mode-coupling i ntegral using the technique 
and approximations discussed in il4l22ll . while adopting the 
mean (zero) value of fci ■ k^, to find 
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The value of each of the integrals for j — 1 — 6 is the 
same when all the \ki\ ~ k. We perform the s-independent 
(case I) trispectrum evaluation for such equal-sided quadrilat- 
eral configurations. Hence, X = = 6 and we 
obtain 



C1234 - ^(^1 + ^2 + ^3 + ^4) X 
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Inserting this into Eq. (HJi for the trispectrum and following the 
approach of |23|, we decompose our delta function as S{ki + 
fc2+fe3+fe4) = / d^K6{ki+k2+K)S{ks+k4-K). Using 
the integral form of the delta functions and the spherical wave 
expansion we perform the integrations over the angular parts 
of (fci, k2, ^3,^4, K), with algebra similar to I.12.,14..241 . to 
give 
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with f. equal to ri for i = 1,2 and r2 for i = 3,4. The ap- 
proximations involved (with respect to angular terms) in the 
ki angular integrals can be made more precise by going to the 
flat-sky limit (elaborated in our detailed paper Il2l[l ). Here the 
K integral gives 6{ri — (7r/2rf ) via the spherical Bessel 
function closure relation. This delta function enables us to 
perform the r2 integral trivially, then ri replaces r2 in the ar- 
guments of j,^ and j,^ . The angular ri and r2 integrals may 
be expressed as [e.g. Eq. 5.9.1 (5) of lEsIl l 
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where we have defined hi^n^ above, along the same 
lines as |23|]. We use the relation [A/S-npo)'^ = 
(2/3)'' {-K/kaf ((n + 3)//fce+^)'' where the Alfven ve- 
locity Va, in the radiation dominated era, is defined as Va — 
So/ (167rpo/3)^/^ « 3.8 x lO^^S^g ill , wifli B„9 = 
(So/IO^^G). From the definition of the rotationally invariant 
angle-averaged trispectrum |26] 
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we separate out the reduced trispectrum ^ (L) (called the 
angular averaged trispectrum in (tE]) from the full trispec- 
trum. We again use the spherical Bessel function closure rela- 
tion to perform the ki integral that yields S{ri—D*) (7r/2r^). 
This facilitates the ri integral that results in ri D* in the 
arguments of j^^ , j^^ and , • The fci, /c2 and integrals con- 
taining a product of a power-law and j"^ can be evaluated in 
terms of Gamma functions (e.g., Eq. 6.574.2 of \2f]). For a 
scale-invariant magnetic index n — >■ —3, we get 
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This gives us the amplitude of the magnetic CMB trispectrum 
sourced by the energy density fi^ of a primordial magnetic 
field. A factor of l/(D*fcG)'*^"^^^ also appears which ap- 
proaches unity for the case n — > —3 of a scale-invariant mag- 
netic field index. We evaluate the magnetic trispectrum for a 
near scale-invariant index n = —2.8, for which this factor is 
^ 1/1500. It turns out that this factor is almost entirely can- 
celled by the increase in value of the k integrals when evalu- 
ated for n = — 2.8 rather than n = —3 112 ill . 

We now compare our magnetic trispectrum with the Sachs- 
Wolfe contribution to the standard CMB trispectrum sourced 



by nonlinear terms in the inflationary perturbations ll23ll28|] . 
More specifically, in the Sachs-Wolfe limit, the dominant term 
of Eq. (64) of Ref. S becomes 



'3 '4 



or. ^ nSW nSW r^sw , 



5.4 X IQ-^'^TNL- 



'hih + iMh + iMk + i)'^' ^^^^ 

Here tjv l and Jn l (below) are standard non-Gaussianity pa- 
rameters and we adopt the standard estimate for the Sachs- 
Wolfe contribution Cf^ 1I23I1 . The factor q which is equal to 
[hih + l)/3(?3 + l)]/[?4(?4 + l)L{L + 1)] is of order unity 
for many configurations. Equation (fT2l i is of the same form 
as Eq. (fTTI) for the magnetic field-induced trispectrum. We 
use the negative-sided limit on tnl derived from searching 
for the CMB trispectrum signal in the WMAP5 data |29], 
tnl > —6000. Magnetic field limits are obtained by tak- 
ing the one-eighth power of the appropriate ratio of trispectra, 
which gives Bq < 16 nG, at a scale of kc = Ih Mpc^^ for 
a magnetic spectral index of n = —2.8. This limit is approx- 
imately 2 times stronger than the Bq < 30 nG upper limit for 
the magnetic energy density bispectrum |12] (taking into ac- 
count the recent estimate of TZ |20]), for the same scale and 
magnetic index. 

We now calculate the trispectrum for the coUinear configu- 
ration [case II]. The full mode-coupling integral 1^-^234 [21] is 
evaluated over all angular terms for the equal-sided collinear 
configuration fci ~ fc2 — ^ ^3 — — ^4. The four-point corre- 
lation of magnetic energy density for the collinear configura- 
tion is found to be 
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Using rt = — 2.8, we compare the collinear configuration four- 
point correlation (, including all angular terms, to ( for case I 
[Eq. |7) that included only s-independent terms. The collinear 
C is similar in magnitude but of positive sign and one then 
expects a trispectrum also of similar magnitude to case I. 

In addition to magnetic energy density, the scalar 
anisotropic stress associated with a primordial magnetic field 
will also act as a separate source for CMB fluctuations - dom- 
inantly in the passive mode IitIi . As we saw in our previ- 
ous work II14II . the magnetic scalar anisotropic stress generates 
^ 10^ times larger contribution to the CMB bispectrum com- 
pared to magnetic energy density. With this motivation and 
using the magnetic trispectrum technique, developed above 
for energy density, we carry out a longer calculation for the 
trispectrum. The temperature anisotropy, sourced via the mag- 
netic Sachs-Wolfe effect by magnetic scalar anisotropic stress 



Hb [defined in Eq. (6) of 111, see also ||I7l|20t]], is 
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where TZp = [-i?^/15] In (Tb/T,.) and Tb and are the 
temperatures at the epochs of magnetic field generation and 
of neutrino decoupling, respectively. 
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For the magnetic scalar anisotropic stress trispectrum, TZ in 
Eq. (4) gets replaced by TZp and C1234 becomes [Ci234]n = 
(nB(fcjnB(fcjnB(fc3)nB(fcJ). fhe full technical details 
of the calculation of the magnetic scalar anisotropic stress 
trispectrum will be presented separately [21]. We give below 
the results considering only the s-independent angular mode- 
couphng terms for equal-sided configurations. In this case 

i.2n+3j.„i.« r(2"/2)(4„ + 3) _ („ + 3)- 
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Here, ^ is a configuration-dependent number that is the sum 
of all s-independent angular terms. This sum involves terms 
like 9ab ~ ka-kb that are constant for a given (fci , fe2 , fca , fc4) 
configuration. Values for ^ range between 2 and 14 for equal- 
sided trispectrum configurations: collinear, square, rhombus 
and tetrahedron. We adopt a typical value ^ ~ 10. This leads 
to a reduced trispectrum 
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We have used Tq ~ 10 GeV (corresponding to the reheat- 
ing temperature) and T^, ~ 10"^ GeV. We see that the ampli- 
tude of the trispectrum sourced by Hb for equal-sided quadri- 
lateral configurations is approximately lO^'' times larger than 
that sourced by fis. Comparison with the trispectrum from 
inflationary perturbations [Eq. [121 gives a magnetic field con- 
sti-aint of 



A/wL < 5 lull from Planck data will imply even tighter 
sub-nG magnetic field upper limits from the scalar anisotropic 
stress trispectrum. Future consideration of magnetic vector 
and tensor modes in the trispectrum is likely to give additional 
constraints on primordial magnetic fields. 

In summary, we have calculated for the first time the CMB 
trispectrum sourced by primordial magnetic fields. The mag- 
netic energy density trispectrum allows us to place stronger 
limits on the primordial magnetic field compared to a simi- 
lar calculation with the magnetic energy density bispectrum 
1 12. [131 . Further, the trispectrum due to magnetic scalar 
anisotropic stress leads to the tightest constraint on large-scale 
magnetic fields of ^ 0.7 nG, approximately 3 times as strong 
as the corresponding bispectrum limit (^ 2.4 nG)[14]. The 
trispectrum's sensitivity is illustrated by the magnetic to infla- 
tionary trispectrum ratio, which is ^ lO'^ compared to 1 for 
the bispectrum (taking f^L ^ 100 and Bq ^ 3 nG). The rela- 
tive contribution of different configurations to the trispectrum 
is different for magnetic compared to inflationary trispectra 
and will be useful to distinguish between them. We also note 
that the magnetic field limit at Mpc scales derived from only 
the scalar magnetic CMB trispectrum is already better than 
the limit (^ 2-6 nG) |7] from the combined scalar, vector and 
tensor modes in the magnetic CMB power spectrum. There- 
fore, the trispectrum turns out to be a new and more powerful 
probe of large-scale primordial magnetic fields. 
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Bn < 1.3 nG, 



(17) 



using the positive-sided limit tnl < 33000 from WMAP5 
data |29]. This is approximately twice as strong as the 2.4 nG 
Bq limit obtained from the Hb bispectrum 1 14] and does not 
assume any particular model of inflation or any relation be- 
tween tnl and /nl- However, for those theories of inflation, 
which lead to tnl — (6/5 fNhf' ifioilsoll . we could perhaps 
use the relatively tighter limits for Jnl- To be conservative 
we take the two-sigma Hmits -10 < /^f ' < 74 on the best 
constrained local Jnl, obtained from searching for the CMB 
bispectrum signal in WMAP7 data 1 3 1 ] . This gives primordial 
magnetic field limits of 



So < 0.7 nG 



and 



Bo < 1.1 nG, 



(18) 



respectively, for the negative and positive /^™' limits. If one 



uses the two-sigma limits for ^^^^ '■^^ ^-^ limit be 

comes 0.6 nG and for 



orthog 
NL 



NL 

it becomes 1.5 nG. However, 



the uncertainties <7f^^ for equilateral and orthogonal config- 
urations are 7 and 5 times larger compared to the local con- 
figuration |31]. Staying with the best determined /j^£°' lim- 
its thus results in sub-nG upper limits on Bq. The expected 
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